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Abstract 


An analysis of the stability of circular rings under a uniformly distributed 
radial load is presented. It was found that, in general, the critical mode shapes 
are a combination of in- and out-of-plane displacements and that they occur at 
loads considerably below the classical (in-plane) critical load. The role of workless 
constraints was also studied. The analysis shows that workless constraints lead to 
an increase in the critical load. 
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The Stability of Circular Rings Under 
a Uniformly Distributed Radial Load 


I. Introduction 

A. Equilibrium Method 

The criterion for stability of the deflected shape of an 
unrestrained circular ring acted on by a uniform distri- 
bution of radial forces was probably first established by 
Levy (Ref. 1) using the equilibrium method. An account 
of this work has been written by Biezeno and Grammel 
(Ref. 2). More recently, Boresi (Ref. 3) sought to refine 
the criteria using the energy method and a deflection 
pattern not bound by the classical assumptions of vanish- 
ing transverse normal and shear strain. In both analyses, 
the form of the disturbance from the original loaded con- 
figuration (the buckling mode) was limited to in-plane 
displacements without twist. 

Most recently, the possibility of out-of-plane buckling 
modes was studied by Wah (Ref. 4). In effect, Wah’s 
analysis was an extension of Timoshenko’s (Ref. 5), which 
had been limited to in-plane disturbances. Wah follows 
Timoshenko’s equilibrium method analysis and intro- 
duces an effective load arising from the ; circumferential 
normal strain into the equations of equilibrium proposed 
by Love (Ref. 6). The resulting equations are uncoupled 
with regard to in- and out-of-plane displacements and 


yield the classical criteria for the in-plane buckling modes 
and a much lower critical load for the out-of-plane modes. 
Wah’s analysis was limited to rings with cross sections 
whose principal axes were oriented parallel to the radial 
and axial coordinate axis (Fig. 1). 

Prior to the work of Wah, a general theory for bending 
and buckling of thin-walled open sections was proposed 
by Cheney (Ref. 7). The method of analysis was essen- 
tially the equilibrium method, in that a set of differential 
equations describing the nonaxisymmetric equilibrium 
the loaded circular state was derived and solved to yield 
the critical load. The cross section was not required to be 
symmetrical as in Wah’s analysis, and the effect of 
warping of transverse cross sections was included by 
incorporating the warping function for straight bars. The 
equations describing the nonaxisymmetric equilibrium 
position were established by requiring that the first 
variation of the change in total potential energy vanish. 

The method followed by Cheney can be shown to be 
equivalent to the energy method in principle. However, 
a number of simplifications were introduced by Cheney 
that might have an appreciable effect if one were to carry 
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Fig. 1. Coordinate system 

out a strict energy method analysis based on Cheney’s 
assumptions for displacements. It is certainly not clear 
that the strain energy density expression used by Cheney 
is as adequate as that expressed in appropriate Lagrangian 
strain components for describing a deviation from the 
loaded circular state. 

B. Energy Method 

The object of the present study is the analysis of the 
stability of a ring acted upon by a constant (in magnitude 
and in direction — see Appendix D), uniformly distributed 
radial load by means of the energy method. This method 
was chosen because it readily allows the nonlinear 
behavior of the structure (prior to buckling) to be 
accounted for. Furthermore, it affords a conceptually su- 
perior framework in which to analyze the axisymmetric 
deviation mode. However, the analysis is not to be 
carried out in the strict sense. It is felt that, for the 
present, it is sufficient to explore the effect of a more 
general energy density expression and the role of the 
axisymmetric mode in a nonlinear analysis while neglect- 
ing the effect of warping. Admittedly, while looking for 
additional terms that might be significant, this one might 
overlook significant warping effects. 

In addition, it is of interest . to obtain the stability 
criteria for circular rings that are not free to deflect 
arbitrarily, but must obey certain kinematic restraints. 
For example, a ring used as a stiffener in a thin- walled 


shell of revolution subjected to external pressure might 
become unstable before the shell becomes unstable. 
The displacements of the ring will depend upon the 
relative stiffness of the shell in and out of the plane of the 
ring, A reasonable method might be to restrict the deflec- 
tions, during buckling of the rings, to those perpendicular 
to the generator of the shell (Fig. 2). Thus, the problem 
to be studied is the stability of a ring acted upon by a 
uniform radial load and subjected to both in- and out-of- 
plane disturbances that may be constrained. However, 
the class of ring cross sections should not be restricted 
to those having principal axes parallel to the radial and 
axial coordinate axis. 

The general stability criteria used here are those 
developed by Langhaar (Ref. 8) and applied by Boresi 
in analyzing the stability of in-plane displacements. 
For a circular ring acted upon by a uniformly distributed 
radial force acting through the line of centroids, there 
will exist, in general, both a uniform radial deflection 
(of the line of centroids) and a rotation of the cross section 
about the line of centroids. If the applied load acting 
over a segment of the ring remains constant in magnitude 
and is either constant in direction or continues to point 
through the original geometric center of the ring when 
the ring is deflected from its equilibrium position, a 
potential can be defined for the applied force and com- 
bined with the strain energy to define a total potential 
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energy. As shown in Ref, 8, the resulting conservative 
system is stable if the second variation of the total 
potential energy is positive definite. Furthermore, the 
system is unstable if the second variation is indefinite, 
negative definite, or negative semidefinite. The critical 
load is defined as that for which the second variation 
becomes positive semidefinite. In particular, Langhaar 
shows that the second variation is positive semidefinite 
if and only if the rank of the associated coefficient matrix 
becomes less than the order of the matrix. Hence, the 
determinant of the coefficient matrix must vanish when 
the second variation changes from positive definite to 
positive semidefinite. This is equivalent to stating that 
one of the stability coefficients (the coefficients of the 
generalized coordinates of the second variation of the total 
potential energy when expressed in canonical form) 
vanishes. 

The energy criterion for stability is essentially that 
suggested by Temple and Bickley (Ref. 9). They argue 
that if a system of disturbances is considered applied to 
the ring, a practical condition for stability is the require- 
ment that the strain energy acquired by the system 
during the disturbance be greater than the work done by 
the primary load (the distributed load) during the dis- 
turbance. The stability limit is reached when the addi- 
tional strain energy is exactly equal to the work done 
by the uniform load. 

In attempting to determine the value of load that would 
cause the second variation of the total strain energy to 
vanish, Langhaar (see Ref. 8) noted that Trefftz recog- 
nized that the existence of a minimum in the second 
variation was equivalent to requiring that the first 
variation of the second variation of the total strain energy 
vanish identically. This requirement leads to a set of 
homogeneous equations whose coefficient matrix is the 
associated coefficient matrix of the second variation. 
Thus, the requirement that a solution exist is that the 
determinant of the coefficient matrix vanish identically 
as before. It can be seen that this form of the stability 
criterion is described by Ziegler (Ref. 10) as the dynamic 
criterion. Clearly, the equation of motion for small dis- 
turbances superposed on an equilibrium configuration 
takes the same form as the equations suggested by Trefftz 
in the limiting case of vanishing frequency. 


of generalized variables that are considered a reasonable 
description of the expected behavior. The form that is 
used is similar to that used for straight beams and has the 
property that the (linearized) transverse normal and 
shear strains vanish at the centroid of the cross section. 
As pointed out in Ref. 8, the resultant stability criteria 
should lead to an upper limit to the true value (compared 
to that for a more general set of displacement variables 
and the same general criteria for stability) of the critical 
load, since the model is stiffer than the actual structure. 

The method of analysis followed in this report ulti- 
mately requires that the expression for the total strain 
energy in the ring be expressible in terms of as few 
displacement variables as possible. However, the number 
of variables must be large enough to give a realistic 
description of the expected behavior. Because it is 
important to retain in-plane displacements with the same 
emphasis as the out-of-plane displacements when at- 
tempting a study of their interaction during buckling, the 
displacements of the ring are described in terms of the 
displacements of the centroidal axis and the (twist) rota- 
tion of the cross section about the centroidal axis. 

In particular, let the displacement components of a 
ring be expressed in terms of generalized variables similar 
to those for straight beams, i.e., transverse cross sections 
translate and rotate as a rigid body. The effect of warping 
is not considered to be important. Thus, with coordinates 
(£, Cj measured with respect to the centroid of the cross 
section (see Fig. 1), a point originally located at 

r = (R + f) e r + £e s 

will move to 

r' = r + [U{e) e r + V(6) e 9 + W(0) ej 

+ [J3(0)e e + K r (0)e r + K z (6)e z ] X (£e r + £e s ) 

As a further condition, let the linearized transverse 
shear strain components vanish at the centroid. In ac- 
cordance with this assumption, K r and K z are determined 
to be 


II. Description of the Displacements 



K z = 


(V-U') 

R 


( )' = 


dO 

de 


As a preliminary step in evaluating the strain energy 
in the ring, the displacements must be expressed in terms 


Thus, the three displacement components u, v, and tv 
are expressible in terms of the four generalized variables 
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U, V, W, and ( ); i.e., 


u = U + £/3 


to — W — /3£ 


R 



(1) 


The strain components consistent with these displace- 
ment forms are derived in Appendix A. 


During the loading process, the uniform radial load p 0 , 
acting through the centroid, is expected to produce a 
uniform radial displacement U 0 and a rotation /?<,. Be- 
cause the loading is axisymmetric, there will be no cir- 
cumferential displacement component. The deviations 
from this uniform (unbuckled) configuration will be 
denoted by U, V, W, and /3; i.e., 


At this stage, the configuration of the ring is described 
completely in terms of the variables U 0 , /?„, A„, B„, C n , 
and D n . 

The strain energy in a ring segment consistent with 
the assumed form of displacement is presented in Ap- 
pendix B. When terms involving deviation components 
in the strain energy expression were expanded, only 
quadratic terms were retained. This restriction applies 
only to the order of the deviation components in any 
given term, because quadratic terms in the deviation 
components that contained U or jB as a multiplicative 
factor have been retained. Furthermore, the circumfer- 
ential strain has been assumed small enough so that 
1 + e e a ~ 1. 

The strain energy expression for a complete ring, 
expressed in terms of the displacement components of 
Eq. (3) and subject to the above limitation, is presented 
in Appendix C. 


u(o) = u 0 + u(e) \ 

m = a . + Mo) > ( 2 ) 

V (tf),W(tf) = V(tf),f(fl) ) 

It will be assumed that the disturbance displacement 
components have the following form: 


(it > p) = X ( An > Bn ) cos n 1 

A / 71 = 0 

-£■= X( c ”) sinn0 

71 = 2 

■ = ^ (D„) cos n 0 


W 

R 




(3) 


The expression for the work performed by the uni- 
formly distributed radial load during a deviation from 
the unbuckled state is derived in Appendix D. It is shown 
that a contribution to the quadratic form for the total 
potential energy results only if the load is assumed to 
be centrally directed. 


III. Analysis of the Unbuckled Configuration 

To apply the stability criteria, the form of the un- 
buckled displacement U 0 and vs p 0 must first be deter- 
mined. Appendices C and D indicate that the total 
potential energy in a complete ring, acted upon by a 
uniformly distributed radial load that maintains its 
original orientation during an arbitrary deviation of dis- 
placement, is given by 


.EAR {(-£)’ (l + ig;)- 


21 rz Ugfio _j I; 


AR 2 R 
\ + G 


AR ! 


PI + 


2\ U„pl 

~E R 


+ 


+ 


E 


n , i n x { i , Ez \ Irz „ , A PoR 

K + 2A o [-R ( 1 + ARt)~ a&Po + TP* ~Ya_ 


2Ro [/»„ (jjfc + f IT + ^ «) - -nr t] + q * ' * + £ (q « “• q ”>} 


(4) 


The generahzed coordinates q 0 = {A 0 , B 0 } T represent the axisymmetric deviation components, whereas q„ = 
{A,„ B„, C„, D„} T represent the out-of-plane deviation components. 
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The equations describing the unbuckled equilibrium 
state result from the requirement that the coefficients of 
the independent deviation components A 0 and B 0 in 
Eq. (4) vanish identically; i.e.. 


/? ( — rr 4- ^ ^0 i ^ ~f~ E „ 2 \ _ Its Up _ 

\AR* + it R + ~k~ J AR 2 R ° ( 5b ) 


Furthermore, the displacement variable p may be elimi- 
nated between Eqs. (5a) and (6b) to give the following 
equation for f(p 0 ): 



( 7 ) 


In terms of the dimensionless variables 


— _U 0 AR 2 x _ n AR 2 _ A + E I 2 r z 

P R I rr ’ ’ Po l n ’ € E AR 2 l rr 

Eqs. (5a) and (5b) can alternatively be written as 



Fig. 3. Rotation of ring vs radial load 


Because Eq. (7) is a cubic equation, it is expected that 
multiple roots will exist for a particular range of the 
parameters e and p 0 (Ref. 11). In fact, as can be seen in 
Fig. 3, a single value of the rotation parameter f exists 
for a given value of load p„ up to a particular value of 
load that depends upon the parameter e. At loads greater 
than this value, the rotation parameter may exist at a 
small positive value f s or at the larger values fi and f 2 
that are approximately equal in magnitude but opposite 
in sign. The corresponding load-displacement diagram 
is presented in Fig. 4. 


The limiting form f or e = 0 and I T z = 0 can most con- 
veniently be seen by examining Eqs. (5a) and (5b). In 
the limit I r * = 0, Eqs. (5a) and (5b) require 



Fig. 4. Load-deflection diagram (v = 1/3, e = 1/400) 
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whenever Eq. (8) is satisfied. 



Fig. 5. Load-deflection diagram for symmetrical 
cross sections iv = 1/3) 


E 1 1-v 
V 2\ 2 V 

These alternative solutions (see Figs. 4 and 5) are 
seen to form a forked branch extending from the point 
p = — E/2\ = —l(v = 1/3) and extending parallel to the 
curves fi, f 2 , and f a corresponding to e^O. The point 
p = —E/ 2A. will be seen to correspond to the limiting 
value of p where multiple values of the rotation param- 
eter f appear as e -» 0. 

To determine the behavior of the solution on the 
parameter e, examine Eq. (6b) near the point p where 
multiple solutions for f(p) appear. According to Burington 
(Ref. 11), the roots of f(p) change character when 



is satisfied. As — p/e ^ 0 for p ^ 0, the roots of Eq. (6b) 
are 


Since e « 1, it is helpful to examine Eqs. (8) and (9) in 
the limit e — > 0. If the three terms in Eq. (6b) are to have 
the same order of magnitude, it is apparent that 

1 + -> 0 as e -» 0 

E 

Hence, when 

1 + -mjL = 8(e) (a 0 + Oi8 + a 2 8 2 + • • •) 

E 

is substituted into Eq. (8), it is apparent that 8 = e v \ and 



Furthermore, the roots of Eq. (6b) corresponding to 
^ = -1 + « 0 € VJ 1 - 2e^ (^J + 3e % + ■ • • J 

( 10 ) 
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are 


where 


U = -2/a 

fc = f. = £(-f )" [l-|e“+(f)V+-.] 

( 11 ) 

For the special case v = 1/3, a 0 — —1.890 and 

p = -1 - 1.890 e V3 (1 + 1.260 + 1.191 + • • •) 

/ 2 (e) = - °-^ 38 (1 + 0.6301 e Vs + 0.3970 e % + • • ■) 


v = p’ + 77ir[' s ( 1+i f)- 2 » + fl 

it follows that V = V(p, /, v, e). 

The function V(p,f ) is determined once the solution of 
Eqs. (5a) and (6b) is obtained. Results of such a calcula- 
tion are presented in Fig. 6 with p cc C7 0 /R as the abscissa, 
and in Fig. 7 with p 0 as the abscissa. The combined 
presentation in Fig. 8 shows only the principal branch 
of the V diagram. In Figs. 6 and 7, it should be noted 
that the value of the strain energy, evaluated at the 
point where the roots of Eq. (5b) change character (as 
defined by Eqs. 10 and 11), is given by 


These relations have been evaluated for e — 1/400 and 
1/1600; the results, together with the value of the load- 
ing parameter at which Eq. (8) is satisfied, are presented 
in Table 1. 


Table 1. Solutions of Eq. (6b) at critical 
point for v — 1 /3 


Parameter e 

Rotation U 

Deflection p 

Load palt'/BIr, 

f = h 

f = h 

1/400 

6.392 

-1.306 

-1.147 

-1.2 77 

1/1600 

9.816 

-1.180 

-1.092 

-1.163 


The load p 0 required at the given value of U 0 /R to 
maintain the configuration f 2 and / 3 must be greater than 
that for /, (see Fig. 4). In particular, the values of the 
loading parameters corresponding to Eqs. (5b) and (8) 
being simultaneously satisfied are 


PoR 3 

El rr 


/=/,, V=Va 


-1 - 1.680 e V3 - 2.646 e % + • • • 


PqR 3 

TTr 


/=/,. V=Va 


-1 - 1.050 6*- 0.265 e % + ••• 


In Eq. (4), the total strain energy should be examined 
as a function of both C/ 0 and p 0 . If Eq. (4) is rewritten in 
terms of the dimensionless p and f variables, and if 


total strain energy = t tEAR 



V 


for the branch f = U and 

-EAR |> - ±al (1 - v)e* + •■■] 

for the branch / = f x . Specifically, 

V(f u e = 1/400) = 1.648 
V(/ 2 ,e = 1/400) = 1.713 



- 3.0 - 2.0 - 1.0 0 


DEFLECTION,-^- 

Fig. 6. Energy vs radial deflection (r = 1/3, € = 1/400) 
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ENERGY, 



Fig. 7. Energy vs radial load (v = 1 /3, € = 1 /400) 



Fig. 8. Principal branch of the strain energy 
vs load-deflection curve 


IV. Axisymmetric Instability of a Complete Ring 

According to Langhaar (Ref. 8), the axisymmetric con- 
figuration specified by U 0 and y3 0 vs p 0 is stable provided 
that the second variation of the total potential energy is 
positive definite; i.e., 

q„ Ko qo + ^ (qj q«) > 0 (12) 

n = 2 

This expression can be put in a more convenient form if 
definitions are given for a complete generalized displace- 
ment vector 


q* = 



and the corresponding coefficient matrix 


“ K„ 0 0 

0 K 2 0 
0 0 K* 


0 “ 
0 
0 


0 0 0 • • ■ K„ 


Hence, the configuration specified by U 0 and f3 0 is stable 
provided that 

q* T K* q* > 0 (13) 


A critical state exists at the load that causes the deter- 
minant of the coefficient matrix K* to vanish. This is 
equivalent to stating that one of the stability coefficients 
(the coefficients of the generalized coordinates when the 
quadratic form is written in canonical form, i.e., the sum 
of squares) vanishes identically. As 

| K* | = | K 0 1 | K 2 1 | K 3 1 • • • | K„ | 

it follows that overall stability can be studied by investi- 
gating the individual stability of each mode. 


This section of the report concerns an investigation of 
the properties of [Ko|. The properties of |K„| will be 
discussed in Section V. 
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In accordance with Appendix C, the determinant of K 0 
will vanish for p and f, satisfying 


1 + 


-rp + , [ 3 ' , -TTr( 1 -¥)’’ = 


(14) 


The coefficients have been simplified with the assumption 
that Izz/AR 2 « 1. Note that if e « 1, the determinant of K 0 
changes from positive definite to negative definite at 

(2 \/E){p) = "I- 


As can be verified (see Ref. 11), Eq. (15) possesses one 
real root, provided that 


T + ir >° 


where 


3q — p 2 

a = -s~ 


7 _ 2p 3 — 9 pq + 27 r 
27 


In the limit e = 0, we must investigate Eq. (14) for the 
two possible configurations that either 

/? 0 = 0, p arbitrary 


When p, q, and r are identified as 


3 (1 — — V) _ 
p 4, (1 + ,) (1 - 2v) 


or 


Pl=~ 


E l rr (, , 2A p \ 
7TTT AR 2 \ ^ E ) ’ 


l + -^£-<0 


As can be easily verified, both choices reduce Eq. (14) to 
1 + 2 \p/E = 0. Hence, a ring of symmetrical cross section 
is unstable with respect to axisymmetric disturbances at 

PoR 3 E_ 

E I rr 2A 

where 


P= ~ 2T 


For e « 1, we expect that the simultaneous solution of 
Eqs. (14) and (6b) would be approximately 


*=--!•■ ' =o(e "> 


The precise values of p and f, satisfying Eqs. (6b) and (14), 
can be obtained as the root of the following cubic equation 
by eliminating p from the set; i.e., 


4f 


E 


) 


-tO-ttf)' 0 < 15 > 


for v = 1/3 


a - . 3 (1 ~ v) 2 _ 3/s 

” n /i i \ /i n \ O/ A 


2 (1 + v) (1 - 2v) 


for v = 1/3 


(1 - v) 2 1 - (1 - v )e 

4ve (1 + v) (1 - 2 V ) 




for v = 1/3 


it follows that, for v = 1/3, 


27 , _ 3 /. e\ 

a 16 ’ 4e 8 ) 


and the criterion is easily satisfied for e « 1. 


In order to determine this single root, let 


1 + ^. = s(b 0 + Sb 1 + S 2 b 2 + ■■■), S = e v 


f = j(fo + hS + f 2 S 2 + •••) 


The values fi can be obtained by substituting the expan- 
sion for f into Eq. (15). A hierarchy of equations results 
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that gives the following values: 


Table 2. Critical deflections for v=l/3 


(1 — v) 2 /4v 

To (1 - 2k) (1 + v) 


U ■ (1 - v)V4v 

ft (1 - 2 v ) (1 + k ) 


(16) 


The corresponding values of b { are most conveniently 
obtained by substituting in the relation 



Critical deflection 

Critical rotation 

Parameter e 




Per 1 1 

Tcrit 

1/400 

-1.303 

7.029 

1/1600 

— 1.181 

10.93 


appear and at the value of the rotation parameter f cor- 
responding to the vertical tangency point on the ft and f 2 
vs p diagram. This interpretation follows from the ob- 
servation that the values of p agree within an error that is 
proportional to e. 


P = f 


1 + ef 



In summary, for a ring with a symmetrical cross section 
e = 0, the deflection states corresponding to that portion 
of the load-deflection curve (see Fig. 5) for which 


obtained by rewriting Eq. (6b). Then, 

+ f [*«-(*-£)(« +£)] + •••} 


so that 


1 - k 3 - 3v - 4 k 2 
fo 4k (1 + k) (1 — 2k) 

(1 - k) 3 3 - 3k - 8k 2 

f 2 0 8k 2 (1 - 2k) 2 (1 + k) 2 


jE / poR 3 

2a ^ \YT^ ’ 


Uo AR 2S 
R I rr , 


<0 


are considered stable with respect to axisymmetric 
disturbances. 


In the general case e ^ 0, the range of load over which 
the ring is stable with respect to axisymmetric disturb- 
ances is slightly larger than | p | ^ E/2\. However, as 
the determinant of K„ changes sign at this critical load, it 
is concluded that the axisymmetric equilibrium configura- 
tion U 0 and R 0 vs p 0 is unstable for loads corresponding 
to (approximately) p = — E/2 a. 


For k = 1/3, we find 

fo = (3/4)* = 0.9086, f t = 1/4, U = f\ = 0.6192 
bo = ~ ~r = - 1.926, b, = - = - 2.271 

*10 Of 0 

so that 

p = -1 - 1.926 8 - 2.271 S 2 + ■ • • 
f = 4- (0.9086 + 0.25 8 + 0.6192 S 2 + ■ • ■) 

O 


V. Out-of-Plane Stability of a Complete Ring 

In this section, the stability of the equilibrium state 
specified by U 0 and R 0 vs p 0 will be discussed with respect 
to nonaxisymmetric deviations. In accordance with 
Eq. (13), the determinant of the overall coefficient 
matrix will vanish if any of the individual 4 X 4 K„ 
matrices possesses a vanishing determinant. 

As is suggested in Appendix C, the analysis of K„ is 
most conveniently described in terms of the elements a,,-, 
where 


These relations have been evaluated for e = 1/400 and 
1/1600; the results are presented in Table 2. 

If the data of Table 1 are compared with those of 
Table 2, it is apparent that the following occurs: the 
axisymmetric deflection state becomes unstable at the 
value of the radial deflection p at which multiple roots 


[ a ij ] 

In forming the determinant of [a j} ], it is important to 
note that all coefficients dij are small (of the order of 
magnitude of e or less) except flu, a ls , and a 33 . Hence, in 
order to work with a well-conditioned matrix, it is nec- 
essary to essentially remove the third column of | |. 
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By appropriate multiplication and subtraction of the 
rows of | da |, it follows that 

| ai > | — a sa | Cij | 


Thus, since e « 1, the elements C i} become 
Cu = e(n 2 - 1)* (l„ + £ + 0(?vs jjj) + 0 (?) 


where the coefficient matrix [C,,] is symmetric and de- 
fined as follows: 


c 

— - /7 — 

/T, 2 

^13 

n 

■ — ft — 

<3_ 

'“'ll 

#11 


022 

#22 


#33 


#33 

C"12 

= ^12 “ 

#23 #13 

j 

£>23 

— #24 

#23 #34 

#33 



#33 

c 

~ &14 

#34 #13 

c 

= #44 

n 2 

W 34 

13 

#33 

^33 

#33 




The significance of removing the third column of 
| da | is apparent if it is assumed that the criterion that 
the quadratic form changes from positive definite to 
positive semidefinite is equivalent to requiring that it 
possess a relative minimum. Hence, the requirement 

8 (qj K » q») = o 


is equivalent to 


K n q„ 0 


C 12 = e(n 2 - 1) I rz + 0(V' 3 1%) + 1„ 0(e 2 ) 

C 13 = n 2 e(n 2 - 1) I„ + 0(e 6 / 3 I* ) + 0(e 2 ) 

C 22 = e (l + ^ + n 2 ^ + 0(eV3 7 % ) + 0(?2) 

C 23 = n 2 e ^1 + -gjPj + 0(e 5/3 Pf) + 0(e 2 ) 

C 33 = n 2 i(n 2 + ^ + + °( ?V3 O + 0(e 2 ) 

If all but the leading terms in the coefficient matrix 
Cjj are omitted, the stability requirement of Eq. (18) 
becomes 



or 

[^ij] Qn 0 

possessing a nontrivial solution. As the first and third of 
this set of four equations is 

A n + nC„ = 0 


rh ( 1 +2^)]=° 

(19) 

For I rz — 0, the critical deflection parameter p C Ht is 
given by 


to order e, it is apparent that the deviation mode shape 
is nearly inextensional in the circumferential direction. 


Since c 33 cannot vanish in the range of interest, the 
stability criterion reduces to 


Cij I — 0 


(18) 


The elements C i} can be expressed simply, since it has 
been noted in Sections III and IV that the range of 
interest of p and f corresponds to 


V = 0(1), f = js 0(1) 


Pent = ~(I - v)n 2 Izz (20) 

and the two additional roots of the quadratic equation 


+ (n 2 - l) 2 


WZ ' 


( 21 ) 


It should be noted that the mode shape corresponding 
to Eq. (20) is an inextensional, in-plane displacement, 
whereas the remaining roots correspond to out-of-plane 
displacements. This uncoupling of the buckling modes 
occurs for I rz = 0 only. 
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For I f g = 7 ^ 0, Eq. (19) is equivalent to the following any curve I z (constant) cannot reach zero without causing 
cubic equation: the parameter T zz — I 2 rs , to vanish. 


+ ;!• [' ' " m + 2 -(”’ + w)]} 

+ (t^t) [<» ! - D’-Jb + <r - “ '»> (1 + w) 

+ 7**-^- (n 2 + 2r)J + (I** — I 2 J (n 2 — l) 2 -^- = 0 

(22) 

It should be noted that the effect of 1 TZ ^0 is reflected 
in the parameter (Igg — 1%) or (1 — l 2 rz /l rr lzz). This is 
also a familiar property of the linear analysis of ring 
structures. Furthermore, Biezeno and Grammel (see 
Ref. 2) indicate that this parameter is always positive. 

The roots of Eq. (22) have been obtained for v = 1/3, 
1/4 ^ I* ^ 3.0, 0 ^ ^ 0.8 I S2 , n = 2, and n = 3. In 

all cases, the least root corresponding to n = 3 has a 
numerical value greater than the least root correspond- 
ing to ft = 2. The remaining roots may overlap. The 
values of the least roots corresponding to n = 2 are pre- 
sented in Fig. 9. It should be noted that the abscissa for 



Fig. 9. Critical deflection vs for an 
unrestrained ring (v = 1/3, n = 2) 


VI. Buckling of a Restrained Ring 

Since the free ring is the exception rather than the 
rule in structural idealizations, it is important to investi- 
gate the effects of restricting the disturbance mode 
shape. In particular, let the radial and axial displace- 
ments of the centroidal axis satisfy the condition that 
there be no component of disturbance displacement 
along an axis inclined at an angle a to the z-axis (see 
Fig. 2); i.e., 

8 • e t = 0 

where 

et = e* cos a + e, sin a 
8-Ue r + fe 2 

In terms of the disturbance parameters, this condition 
requires that 

A n sin a + D n cos a = 0 

Because D n is now dependent upon A n , the criterion 
for stability derived in the previous section must be 
modified. If the analysis is restricted to forces that main- 
tain their original line of action after buckling, the only 
modification that is necessary to the quadratic form is 
the elimination of D n as a variable. In this way, the 
holonomic character of the system is retained. 

With regard to the nth mode, the substitution of 

D n — —A n tan a 

changes the quadratic form in the total potential energy 
expression of Eq. (4) to 

(A n , B n , C„) K£* ) | B* | 

where K£ 0) = [k^f ], and 

= fcu — 2 k 14 tan a + k 44 tan 2 a 

KV = *22 KV = *33 

= k r2 — k 2i tana = 
k[f = k 13 - k 3i tan a = 

*23 *23 * 32 
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The variable C„ is essentially eliminated by removing 
the third column. In terms of the coefficient C i} defined 
in Eq. (17), the determinant of can be shown to 
vanish for 

|%|=° (23) 

where 

Bn = Cn — 2 C a3 tana + C 33 tan 2 a \ 

Bi 2 = C 12 — C 23 tana = B 21 > (24) 

B 22 1 C 22 / 

For a = 0 and I r * 0, the criterion for stability is that 
p satisfy 

C u c 2 , - C 2 2 = 0 (25a) 


Taking only the leading terms in the expression for the 
elements C^, we obtain the following quadratic equation 
for p: 



+ 



^1 + n 2 


GJ 

Wr 


+ 2 vn 2 Ig g 


) 


+ (La — l 2 „) + n 2 I z!! 


GJ 

Elrr 


= 0 


(25b) 


In the general case a^O and I rs ^0, Eq. (23) re- 
quires that 

(GnC 2 2 c\ 2 ) 2 tan a (C 13 C 22 G 12 C 2S ) 

+ (C 22 c 33 - C 2 J tan 2 a = 0 (26a) 


Taking only the leading terms in the expressions for the elements C ih we obtain the following quadratic equation 
for p: 


2v 

n 2 




GJ 


4vn 4 I r 


^Et r +2vn%1! ~ n 2 - 


■ tana 


+ 


+ - 1% + n 2 (Lz - 21 rz tan a + tan 2 a) j = 0 


(26b) 


Note that the above analysis presupposes that tan a 
= 0(1), since A n is used as a variable. For restraint con- 
ditions corresponding to a at approximately 90 deg, D„ 
would be a more convenient variable, and the analysis 
would proceed in a manner analogous to that above. 

The roots of Eq. (26b) have been obtained for v = 1/3, 
1/4^ 3.0, -0.87*5^ Irz ^ 0.8T;, n = 2, and a = 0, 

15, 30, 45, and 60 deg. The least of the two roots is 
shown plotted in Figs. 10-14. A summary of the results 
for a symmetrical, restrained ring is presented in Fig. 15. 
Note that the effect of restraint for unsymmetrical cross 
sections lm¥= 0 is a general increase of the critical de- 
flection parameter. 


VII. Discussion 

To interpret the significance of the various stability 
criteria established in Section VI, the load-deflection 


diagrams should be reviewed (see Figs. 3, 4, and 5). For 
a ring with a symmetrical cross section, the radial deflec- 
tion C7 0 is directly proportional to the radial load p 0 , and 
the deflection increases with no rotation of the cross 
section /3 0 = 0. At the load corresponding to a radial 
deflection p = —E/2\, a multiple-valued deflection con- 
figuration appears with no jump in the rotation. For 
values of load corresponding to a radial deflection greater 
than — E/2\, three different rotations exist for a given 
radial deflection, namely, 


Bo = 0, 



~E Irr 

~kTT!! ar 2 



As shown in Section IV, the point at which this multiple- 
valued configuration appears is unstable with respect to 
axisymmetric disturbances. Hence, it is unlikely that a 
ring of symmetrical cross section could be loaded above 
this point. 
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DIMENSIONLESS AREA MOMENT, 



-3.0 -2.0 -J.0 0 

CRITICAL DEFLECTION PARAMETER, p 

' r cnf 


Fig. 10. Critical deflections for a restrained 
ring (a = 0 deg, n = 2) 



-4.0 -3.0 -2.0 -1.0 0 

CRITICAL DEFLECTION PARAMETER, p .. 

r cnf 

Fig. 1 1 . Critical deflections for a restrained 
ring (a = 15 deg, n = 2) 

For a ring with an unsymmetrical cross section, the 
radial deflection U 0 is always accompanied by a rotation 
/3 0 for all load values. If the ring were loaded by a con- 
strained, radial deflection device, the radial deflection 


and the rotation would be single-valued functions of 
the load up to a load value corresponding to a radial 
deflection slightly greater than p = — E/2X. At this point, 
three values of rotation can exist for the same radial 
deflection. Hence, there may be a jump in the rotation 
and a consequent increase in the radial load required to 
hold the ring in the new position. The value of the rota- 
tion that appears at this point is opposite in sign to that 
which accompanied the radial deflection from the un- 
loaded state. As the radial deflection increases from this 
point, two values of the rotation parameter increase in 
magnitude with the load (the f t and f z branch), whereas 
the third value tends to decrease in magnitude with 
the load. 



Fig. 12. Critical deflections for a restrained 
ring (a = 30 deg, n = 2) 


As shown in Section IV, the determinant of K 0 changes 
from positive to negative definite at a load correspond- 
ing to [ p [ slightly in excess of E/2\. Hence, one must 
conclude that an axisymmetric deflected state cannot 
exist for loads in excess of approximately — E/2\. The 
classical critical load is, therefore, impossible to achieve. 

Strictly speaking, the actual value of the radial load 
corresponding to the load parameter p 0 must be found 
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- 3.0 - 2.0 - 1.0 0 

CRITICAL DEFLECTION PARAMETER, p 

r cnt 


Fig. 13. Critical deflections for a restrained 
ring (a = 45 deg, n = 2) 

by considering the length of the centroidal axis in the 
deformed configuration. The form of the energy principle 
used defines the distributed force per unit length of 
centroidal axis in the undeformed state as p 0 and the 
correspondence rule indicates that p 0 is also the force 
per unit length in the deformed configuration. Hence, 
the total load acting on the ring in the deformed state 
is equal to the total load acting on the ring in the 
undeformed state multiplied by the ratio of the unde- 
formed periphery to the deformed periphery. Consider- 
ing that the linearized circumferential normal strain e 0 o 
was neglected in comparison with unity in computing 
the total strain energy, there is essentially no distinction 
between either the periphery or the total load in either 
the deformed or the undeformed state. 

The results of the stability calculations made in Sec- 
tion V apply only to very thin rings. In particular, 
Eqs. 19-22 were written with the assumption that only 
the leading terms in the coefficients C l} need be retained. 
This implies, in general, that terms of order e.' A = ( l„ / 
AR 2 ) Va were negligible in comparison with unity for 
T r z = Irz/lrr = 0(1). For I r z = 0, the assumption is that 
O(e) is negligible in comparison with unity. Hence, if 



- 3.0 - 2.0 - 1.0 0 

CRITICAL DEFLECTION PARAMETER, p 

crit 


Fig. 14. Critical deflections for a restrained 
ring (a = 60 deg, n = 2) 

improved accuracy is warranted for l rz ^= 0, the tech- 
nique employed in Section V must be replaced by a 
search technique wherein successive values of p and the 
corresponding values of f x would be substituted into 
| Cij | until ] C ,- j | vanished. 

This procedure would require the simultaneous solu- 
tion of Eq. (6b) (or a table of values of f vs p) for each 
step. It is only for the limiting case described in Sec- 
tion V for e « 1 that the terms in ] C i} | do not involve 
f L ; hence, simple cubic and quadratic equations, Eqs. (21) 
and (22), result. An analogous procedure would also have 
to be undertaken for the constrained ring, as the same 
approximations were made in Section VI as in Section V 
regarding e « 1. 

Furthermore, the radial load parameter p, in the 
discussion above, was considered to be equal to the load- 
ing parameter p 0 R 3 /EI rr of Eq. (6a). Since the value of f u 
corresponding to the critical value of p, has been shown 
to be near 0(1) /e A (see Section III), it follows from 
Eq. (6a) that taking p « p 0 R 3 /EI rr also involves regard- 
ing e I/3 as negligible in comparison with unity. Hence, as 
in the evaluation of I C,-,- 1, the explicit evaluation of 
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CRITICAL DEFLECTION PARAMETER 



0 30 60 


SLOPE OF RESTRAINT a, deg 


/ vs p is only required when the assumption e Va « 1 is not 
warranted. 


The results for the nonaxisymmetric modes based 
upon Eq. (19) indicate that this analysis does not recover 
the classical result achieved by Timoshenko (see Ref. 5). 
The classical result is interpreted here to mean the 
critical load corresponding to in-plane displacements 
only. To obtain such a mode, one must examine the 
limiting case I n = 0. Following Eq. (20), and noting 
that p = p 0 R?/El„ for I r * = 0, and fi„ = 0, we see that 


M!\ 4(1 - vY 

Eh,) crit (l + v)(l-2v) 

= 4.00 (r = 1/3) 


n = 2, v arbitrary 


= 3.60 (v = 1/4) 


The minimum value of this expression is 3.56 and occurs 
for v = 1/5. This formula is essentially the result of set- 
ting C u = 0. 


When comparing the minimum value with the classical 
value (3.00), one should note that the result reported by 
Boresi, as well as Wempner and Kesti (see Ref. 12), was 
that the classical value corresponds to a hydrostatic 
type of loading in which the load at a point on the ring 
remains normal to the centroidal axis during buckling. 
The type of loading studied here (see Appendix D) is 
that identified by Wempner and Kesti as constant — the 
orientation of the loading does not change during buck- 
ling. The value of the critical load for constant loading 
was found to be 4.00 by Wempner and Kesti using the 
equilibrium method. Note that the identical result is 
obtained here for v — 1/3, and that Wah obtained the 
value 3.00. 

The results of calculations based upon Eq. (21) for 
nonaxisymmetric buckling of symmetrical sections might 
be qualitatively compared with the results of Cheney 
and Wah. The intention here is that the dependence of 
the out-of-plane buckling loads on the postbuckling be- 
havior of the loading be noted rather than that any 
quantitative comparison be made. As can be seen in 
Table 3, the results given for the constant load are con- 
sistently about half the values obtained by Cheney and 
Wah for the hydrostatic type of loading. 


Fig. 15. Summary of critical deflections for a 
symmetrical, restrained ring 


Salient characteristics of the numerical results of this 
analysis are presented in the concluding statements that 
follow. 
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Table 3. Comparison of critical load for 
symmetrical sections ((» = 0) 


Area moment 

/»* 

Results 

Present 

<*=l/3) 

Ref. 4 

0.25 

-0.5948 

-1.467 

0.50 

-0.6800 

-1.558 

1.00 

-0.8258 

-1.688 

2.00 

-1.043 

-1.842 

3.00 

-1.194 

-1.929 


The decreasing slope of the curves I z (constant) should 
not necessarily be interpreted as a destabilizing effect 
due to I rz 0 (see Fig. 9). A more instructive cross plot 
might be made for a given cross section showing the 


variation of the critical load with inclination of the 
principal axis. In fact, such a cross plot would probably 
be similarly inclined, but not as sharply because the 
parameter I* changes with inclination. The same is true 
for the restrained ring. With the latter, however, the lack 
of symmetry of the results with respect to the parameter 
I r e indicates that there is a preferred direction for a given 
cross section having a greater critical load. 

Finally, it should be noted that the effect of restraint 
(see Fig. 15) for symmetrical cross sections is true in 
general; that is, the critical load for the restrained ring 
is greater than the critical load for the corresponding 
unrestrained ring. Thus, the behavior of curves for I rz 
(constant) is qualitatively similar to that of the curves of 
Fig. 15 for all values of the parameters studied in 
Figs. 10-14. 
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Appendix A 

Strain Components in Circular Polar Coordinates 


Let dR = e r dr + e« rd e + e* dz be a differential line 
segment in an elastic material that becomes dR' when 
the material is deformed with the displacement field 

8 (r,0,z) = u(r,6,z) e r + v( r,6,z)e 0 + w(r,0,z)e z 

As 

dR' = dR + d5 


2E r * = W* + 10 r + (u r Ug + V r Vg + IVr iOg) 

2 Egg — Vg + — Wo + — [Ug (u 0 — v) + Vg ( U + Vo) + Wo Wg] 

r r 

where l, m, and n are the direction cosines of the line 
segment dR, i.e., dR = j dR |(Ze r + me, + ne z ). 


= e r [(l + u r )dr + (tt s — v)d6 + u z dz\ 

+ e 8 j^o r dr + rdd ^1 H — ^ + Vg d z J 
+ e z [w r dr + Wo d0 + (1 + w z )dz ] 


In terms of the generalized coordinates U, V, W, and 
fi, the strain components become 



and defining the components of the Lagrangian (Green’s) 
strain tensor as 

= l * E rr + m*Eoo + n?Egg + 2 Im E re 

2(dRy 

+ 2 InErz + 2mnE„ 

it follows that 


E„ = (l + f) + -V) + ({/>■ + f)J 

+ (w - zpy | 

re.. = u + Ve = V + 03 + -£-V'- jj-U" - j^W" 


E„ = U r + (it* + Vl + Wl) 



Eoo = “ (« + Ve) + -jp- [(Ue ~ v) 2 + (u + V„) 2 + W\] 
Egg = Wg + -j- (u 2 e + v 2 z + wl) 

2 E r0 = v r + -jj- (Ue — v) 

+ -i- [u T (u e — v) + Vr (tl + Vo) + W r tt> s ] 


2Ero =-f (V + W) + eoo 

W 

2E r g=-^ r (V-U') 

1 / W'\ U'—V W' 

2Eog = - T \B' + - w Ht-pO + P — r e«-R- 
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Appendix B 

The Strain Energy in a Ring Segment 


The strain energy # sC 0 2 in a ring segment can be 
written, following Fung (Ref. 13), as the integral of the 
strain energy density p 0 W over the volume of the ring in 
its undeformed configuration, i.e., 

[ p„W(E i3 )dV„ 

Jv„ 

where 

_ d(poW) 

ij dE i} 

are the components of the Kirchhoff stress tensor. For a 
linearly elastic, isotropic material, 


Sij — \E/ac Sij + 2 GEij 

where \ and G are the Lame constants of the material, it 
then follows that the strain energy density takes the form 

p 0 W(E ti ) 

+ A ( ErrEoo + EzzEeo + E rr E^z) 

+ G(El 0 + + El + El + El + El) 

where E = A + 2G. 


In terms of the generalized coordinates U, V, W, and /?, the individual terms in the strain energy expression become 




W" \ 
R ) 


h 


^( U iJ 




R h 


±iin 

R Lz 


-) 


+ 2 




+ 


■^+D’d 


TT I* 

ii. + a - 

R P Rlrr 


) 


JrdAE" (E rr + E„) = ARp j^— + V ' + ~ ^ 

L ( ff + IV + i / Ziyi , ar u + v' r/wy 

2 AR 2 y R )' 2 \ R ) 2 R [\ R / 

(V-,).] 


+ ■ 

+ 


f ;i AE„ E„ - lAR,. |> + (f)' + (2=*)' 
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4 ^ (E , +Ey= x^ + -)% A „[(-)V(v i -y] 

x[(4-) ! ( i+ ^) + K 1+ w)-^] 

where terms in U', V, and W' of order higher than quadratic have been discarded, and e „ 0 and [i 2 have been neglected 
in comparison with unity. Note that higher-order terms have been retained if they contain U or /? as a factor. The fol- 
lowing area properties have been used: 
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Appendix C 

Strain Energy Relations for a Complete Ring 


If the form for the displacement components given in Eqs. (2) and (3) are substituted into the individual strain 
energy contributions of Appendix B, the following is obtained:- 


r*9 c 
J 0=0 Ja 


r de dA El + E* ) = ttAR/31 <J /? 2 + 4 R 0 B 0 + 6B 2 + 3B 2 + - J- D 2 + -j- (C n + nA 

n - 2 


-n 


n) 2 | 


r dO dA E rr E~ 


r 2 7r r 

Je = o Ja 


r dO dA e 2 = 


2 ,AH [(-£■ + A.)‘(l + -jjfr) - + A.) (A + ft) + 4gr (A + ft)’] 

+ ttAR JT (A„ + nC„) 2 + -^fr (n 2 - l) 2 A » + -^ (B» + " 2 Dn) 2 + ^ (" 2 - 1) A n (B„ + n 2 D„) 1 

2 

£ {/ MdAE i. = £L ,MdAe '“ + * ar -t 2[< c - + nA ->* + 

n-2 

- 7§F "< B ” + D ’> < c - + » A -> (x “ T5r) + W"’ < B * + D -) ! (x- + 

n r 217 

r d0 dA E $$ (E rr + E, z ) = irAR)A~(i3 0 + B 0 ) 2 + 2 A 0 (/? 2 + 2/? 0 B„) 

v 

+ [ n2 Dl + (^» + wA„) 2 + 2B 2 ] + 2/3 0 ^ B„(A„ + nC„) 


+ 


T S [ ^ + C ” )2 + W (B " + Dn) 2 + n2D nl} 


Z*2?r /• 

/ rdOdA(El + E*J = 

Jo = 0 7 A 

S iT? < B * + D -) ! + ir (ft + D,.) [ l„ (C, + »A.) ^ h±hhi - p, I zJJhht ) 

- nD * (1- - * - Hf" 1 " )] + ^ ” !D * + < c * + ,,A -) 2 ] 

x [(!)’ (! + w) + « (l + ^r) - 
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Hence, the strain energy in a complete ring can be written as a sum of contributions from each order n of the devia- 
tion displacements. For reference, it is convenient to write the strain energy in the following form: 


Strain energy ttEAR^ r ^ ^1 + AR ^j AR2 R + AR2 PI 

_ 2 a_ Uofi. A + G 04 i o a f 02 _l Hz / 1 _L ill \ _ isi o "1 

+ E R + E + 2A ° L E + R V 1 + AR 2 ) AR 2 J 

■ or fo / hz I 2a U n A + E \ _ l rz U 0 

+ 2B °(_^0\AR 2 + E R + E P°) AR 2 R J 

+ q 0 T K„ q 0 + ^ (q* K n q„) j 

n-2 ' 


The vectors q 0 and q„ and the matrix K 0 are defined as follows: 



The 4X4 matrix K n is most conveniently represented as 

K|i [^i;] ) 2kij Q-ij 


where the elements a i; - are 

flu = 1 + € £(» 2 - l) 2 f^ + P n ' 2 A + -jr p 2 n 2 e 2 + 2n 2 e 2 f^ ^ G - f - Jpf^j 

a 12 = en 2 l ra + (P ~ f) + ? f % (P h ~ f % J j 

<*13 = »[ 1 + + -f e 2 p 2 + 2e 3 -| Z 2 z /p] 

a u - n 2 I„e n 2 - 1 - e(p - f ) - A - g G + e 2 / -jr (p Z* - fPJ J 


-1", 2A p , „ GJ 2 /-/* ,JT 

a 22 — e 1 + -p- + n 2 — h e n 2 1 p j — h f I 


J * -I_ T* N 

T I ‘rrz A rrr 


_ j_ !L g /2 72 



e±+^(p-f) + e> 
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where 


and 


d2i 


_ , G] , _/-T A + G , _ /* , , r It rr -R'i% , G fF , \ 

= n 2 e 7 1 + ttt help L — ^ bp- b fl ,* — 1- ^t/I 2 * j 


EL 


\ 


E 


RL 


d-SR tl 2 + 
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Appendix D 

The Work Expression 


The work performed by a uniformly distributed radial 
load during a deviation from the unbuckled state depends 
upon the orientation of the forces acting on a line segment 
during the displacement. The most realistic distinction 
seems to be that between (1) a force that retains both 
its magnitude and direction, and (2) a force that is 
constant in magnitude but whose line of action passes 
through the undeformed center of the ring. 

In general, a differential force dF acting on a line 
segment R dd of the line of centroids will perform work 
dW e during a displacement r 2 — r, (Fig. D-l) where 

dW e = f" dF-dr 

Jr — pi 

For a constant load, dF x = p 0 Rd8 e' 0) 



Fig. D-l. Sign convention for work expression 


For a centrally directed load, dF 2 = p 0 Rd$e r 

If dr = e r dr + e» rdd + e« dz, the differential work 
expressions become 

dW = p 0 RdB e< 0) (r 2 — r t ) 

dW ( e 2) — p 0 Rde (r 2 — r ± ) 

Furthermore, if r 2 — 1 ^ = Ue l r 0) + Ve< 0) + We™, it 
follows that 

dW<» = p 0 Rd8U 
Also, with r\ = (r t + U) 2 + V 2 + W 2 , 
then 

r 2 -r 1 = U + ~(V 2 + W 2 ) + ■■■ 

LTi 

and 

dwy = p 0 Rdo^u + ~^r(V 2 + W 2 ) + — J 

Finally, using the deviation displacement forms of 
Eq. (3), we find that the total work expressions become 

= 2nr p 0 R 2 A„ 

Wf > = %r Po R*A 0 + (Cl + D% 
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Nomenclature 


A area of cross section 
An Fourier coefficient of radial deflection 
an two times the elements of the stiffness matrix 
Bij elements of a coefficient matrix 
B n Fourier coefficient of tangential rotation 
Cij elements of a coefficient matrix 
C n Fourier coefficient of tangential deflection 
D n Fourier coefficient of axial deflection 
E Young’s modulus of material 

E Lame constant of material plus twice the shear 
modulus of material 

f nondimensional rotation 
G shear modulus of material 
Iij area moment 
1 dimensionless area moment 
J area polar moment 
K„ stiffness matrices 
kij elements of the stiffness matrices 
n Fourier wave number 
p 0 intensity of distributed radial load 


p nondimensional radial deflection 
q„ generalized displacement vector 
v radial coordinate 

R radial coordinate to centroid of section 
U radial component of displacement 
V tangential component of displacement 
W axial component of displacement 
z axial coordinate 
a constraint parameter 
fi tangential rotation component 
e dimensionless parameter 
e slenderness ratio 
£ centroidal axial coordinate 
6 tangential coordinate 
K r radial rotation component 
k z axial rotation component 
X Lame constant of material 
v Poisson’s ratio 
£ centroidal radial coordinate 
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